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Abstract 



1 Introduction 



We have examined the dust photoelectric heat- 
ing in the intergalactic medium (IGM). The 
heating rate in a typical radiation field of the 
IGM is represented by r pc = 1.2 x 1CT 34 
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K) -1 / 6 , where T> is the dust-to-gas mass 
ratio, riH is the hydrogen number density, Jl is the mean 
intensity at the hydrogen Lyman limit of the background 
radiation, and T is the gas temperature, if we assume the 
new X-ray photoelectric yield model by Weingartner et 
al. (2006) and the dust size distribution in the Milky Way 
by Mathis, Rumpl, & Nordsieck (1977). This heating rate 
dominates the HI and Hell photoionization heating rates 
when the hydrogen number density is less than ~ 10~ 6 
cm~ 3 if V = 10~ 4 which is 1% of that in the Milky Way, 
although the heating rate is a factor of 2-4 smaller than 
that with the old yield model by Weingartner & Draine 
(2001). The grain size distribution is very important. 
If only large (> 0.1 /im) grains exist in the IGM, the 
heating rate is reduced by a factor of ~ 5. Since the 
dust heating is more efficient in a lower density medium 
relative to the photoionization heating, it may cause an 
inverted temperature-density relation in the low density 
IGM suggested by Bolton et al. (2008). Finally, we have 
found that the dust heating is not very important in the 
mean IGM before the cosmic reionization. 
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Dust grains are formed at the end of the stellar life; in the 
stellar wind of asymptotic giant branch stars (e.g., Fer- 
rarotti & Gail 2006), in the stellar ejecta of supernovae 
(e.g., Nozawa et al. 2003, Rho et al. 2008), and so on. 
Some of the grains grow in molecular clouds (e.g., Draine 
1990). Some of them are destroyed by the interstellar 
shock (e.g., Williams, et al. 2006). Some of them may 
go out from the parent galaxy and reach the intergalactic 
medium (IGM) (e.g., Aguirre et al. 2001a,b). 

The IGM is the medium between galaxies, and it oc- 
cupies almost the whole volume of the Universe. The 
mean density of the IGM is as low as 10~ 7 -10 -4 cm~ 3 . 
As found by Gunn & Peterson (1965), the IGM is highly 
ionized after the cosmic reionization epoch (the redshift 
z ~ 6-10; Loeb & Barkana 2001, Fan, Carilli, & Keat- 
ing 2006). Thus, its temperature is ~ 10 4 K. The IGM is 
filled with the ionizing ultra-violet (UV) and X-ray back- 
ground radiation which is produced by QSOs and galaxies 
(e.g., Haardt & Madau 1996). 

A significant amount of metals is found in the IGM 
(e.g., Aguirre et al. 2001c). Multiple supernova explo- 
sions (SN) caused by an active star-formation in galaxies 
can eject the metal elements to the IGM. However, Fer- 
rara, Pettini, & Shchekinov (2000) showed that the metal 
enrichment of the IGM by SN explosions is limited in rel- 
atively small regions around star-forming galaxies and an 
additional physical mechanism is required to explain the 
observed global enrichment of metals in the IGM. Dust 
grains expelled from galaxies by the radiation pressure 
due to stellar light and by the galactic wind due to multiple 
SNe may contribute to the metal enrichment in the IGM 
(e.g., Aguirre et al. 2001a,b). Bianchi & Ferrara (2005) 
showed that relatively large (> 0.1 /im) dust grains are 
not completely destroyed and reach a significant distance 
(a few x 100 kpc) although the amount of this intergalac- 
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tic dust is too small to make a detectable extinction. 

Infrared (IR) emission from dust grains in the IGM 
surrounding edge-on galaxies has been already detected 
(e.g., Alton, Davies, & Bianchi 1999, Bendo et al. 2006). 
Moreover, IR emission from dust in the IGM accumu- 
lated from the distant Universe may affect the cosmic 
far-infrared background and the cosmic microwave back- 
ground (Aguirre & Haiman 2000). Emission signature 
from dust even at the epoch of the cosmic reioniza- 
tion may be detectable with a future satellite observing 
the cosmic microwave background (Elfgren, Desert, & 
Guiderdoni 2007). 

Xilouris et al. (2006) found a significant reddening of 
galaxies behind the M8 1 group IGM which is detected by 
HI 21 cm emission (e.g., Yun, Ho, & Lo 1994). Their 
measurements imply the dust-to-gas ratio in the M81 
group IGM is a factor of 5 larger than that in the Milky 
Way. Such a large amount of dust in the IGM may be 
ejected from M82 by its intense starburst activity (Alton 
et al. 1999). 

Dust in the IGM affects results from the precision cos- 
mology. Indeed, high redshift SNe la are dimmed by dust 
in the IGM, and then, the observational estimate of the 
distance to them and cosmological parameters become 
ambiguous (Goobar, Bergstrom, & Mortsell 2002). Fur- 
thermore, future investigations of the 'equation of state' 
of the Dark energy will be affected by the extinction of 
the intergalactic dust even if its amount is too small to 
affect the conclusion of the presence of the Dark energy 
(Carasaniti 2006, Zhang & Corasaniti 2007). 

Dust in the IGM also affects the thermal history of 
the IGM. In the intracluster medium, dust grains work 
as a coolant because they emit energy obtained from gas 
particles collisionally as the thermal IR radiation (Mon- 
tier & Giard 2004). Such emission from some nearby 
galaxy clusters can be detectable with the current and 
future satellites for the IR observations (Yamada & Ki- 
tayama 2005). Dust grains in the IGM also work as a 
heating source via the photoelectric effect (Nath, Sethi & 
Shchekinov 1999). Inoue & Kamaya (2003, 2004) pro- 
posed a possibility to obtain an upper limit of the amount 
of the intergalactic dust based on the thermal history of 
the IGM with the dust photoelectric heating. 

In this paper, we revisit the effect of the dust pho- 
toelectric heating in the IGM. Recently, Weingartner et 
al. (2006) revised the model of the photoelectric yield 



of dust grains. They included a few new physical pro- 
cesses; the photon and electron transfer in a grain, the 
photoelectron emission from the inner shells of the con- 
stituent atoms of grains, the secondary electron emission, 
and the Auger electron emission. These new features re- 
duce the photoelectric yield for moderate energy photons 
of ~ 100 eV but enhance the yield for high energy pho- 
tons of > 1 keV. In particular, we explore the effect of 
the new yield model on the photoelectric heating by the 
intergalactic dust in this paper. 

The rest of this paper consists of four sections; in §2, 
we describe the model of the photoelectric effect. In §3, 
we compare heating rates of the photoelectric effect with 
those of the photoionization in the IGM. In §4, we dis- 
cuss the implications of the results of §3. Final section is 
devoted for our conclusions. 

2 Dust photoelectric effect 

2.1 Grain charging processes 

To examine the photoelectric effect, we must specify the 
charge of grains which is given by the following equation 
(Spitzer 1941, Draine & Salpeter 1979): 

— Rj + 7 (i) 

i 

where is the grain charge in the electron charge unit, 
Ri is the collisional charging rate by i-th charged particle 
(hereafter the subscript "i" means "z-th charged particle"), 
and i? pc is the photoelectric charging rate. We consider 
only protons and electrons as the charged particle. 

2.1.1 Collisional charging rate 

The collisional charging rate by i-th charged particle, R i7 
is expressed as (e.g., Draine & Sutin 1987) 

R i = Z i s i n i l ai(a,Z d ,Zi,Vi)vif(vi)dvi, (2) 
Jo 

where Zi is the charge in the electron charge unit, s, is the 
sticking coefficient, rij is the number density, Vi is the ve- 
locity, G{ is the collisional cross section depending on the 
grain radius, a, both charges, and the velocity, and f(vi) 
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is the velocity distribution function. If the grain and the 
charged particle have the charges of the same sign, the ki- 
netic energy of the particle must exceed the grain electric 
potential for the collision. Otherwise, the collisional cross 
section is zero. We simply assume Sj is always unity. 

Now, we introduce the dimensionless cross section, 
5» = cfi/ira 2 . If we neglect the "image potential" re- 
sulting from the polarization of the grain induced by 
the Coulomb field of the approaching charged particle 
(Draine & Sutin 1987) and we assume the Maxwellian ve- 
locity distribution for the particle and the spherical grains, 
we obtain 

/ aiVif(vi)dvi = g(x), (3) 

Jo V * m i / 

and 

_ f 1 - x for Z d Z, < 
9[x) ~ \ exp(-ar) for Z d Z t > ' (V 

where k& is the Boltzmann's constant, T is the gas tem- 
perature, m, is the particle mass, and x — e 2 Z d Zi/akBT 
(Spitzer 1941). 

In fact, the "image potential" works to enhance the col- 
lisional cross section (Draine & Sutin 1987). Although 
the effect becomes the most important for grains with 
an around neutral charge, it quickly declines for highly 
charged grains. Indeed, for the charge ratio of Z d /Zi < 
—3, which is satisfied in our case as found below, the in- 
crement factor for the cross section by the effect of the 
"image potential" is less than 1.5 (Draine & Sutin 1987). 
Therefore, we neglect the "image potential" in this paper. 

2.1.2 Photoelectric charging rate 

The photoelectric charging rate is given by (e.g., Draine 
1978) 

r°° At:J v 
R pc = ira 2 Q v {a)Y v {a,Z A )——^-dv, (5) 
Jo < u ' 

where Q v is the absorption coefficient of grains at the fre- 
quency v, Y v is the photoelectric yield, J v is the mean 
intensity of the incident radiation, and h is the Plank con- 
stant. For Q u , we adopt the values of "graphite" and "UV 
smoothed astronomical silicate" by Draine (2003). If the 
photon energy is smaller than the threshold energy of the 



photoelectric emission, e.g., the ionization potential or the 
work function, the yield Y v = 0. 

We adopt a sophisticated model of the photoelectric 
yield by Weingartner & Draine (2001) and Weingartner 
et al. (2006) in this paper. The model of Weingartner 
& Draine (2001) (hereafter WD01 model) takes into ac- 
count the primary photoelectron emission from the band 
structure of grains, a small-size particle effect, and the en- 
ergy distribution of the photoelectron. On the other hand, 
Weingartner et al. (2006) (hereafter W+06 model) add the 
primary photoelectron emission from inner shells of the 
constituent atoms of grains, the Auger electron emission, 
and the secondary electron emission produced by primary 
electrons and Auger electrons. The transfer of photons 
absorbed and electrons emitted in a grain is also taken 
into account. The detailed procedure of the yield calcu- 
lations is referred to the original papers of Weingartner 
& Draine (2001) and Weingartner et al. (2006). Fig. 1 
shows comparisons between the WD01 and W+06 mod- 
els. The reduction of the W+06 yield around 100 eV is 
due to the effect of the photon/electron transfer in a grain. 
The W+06 yield exceeds unity for some cases because of 
the Auger and secondary electrons. 

We have to note that there is still a large uncertainty of 
photoelectric yield models because of our insufficient un- 
derstandings of the nature of small-size particle effect as 
well as the lack of experiments. Abbas et al. (2006) re- 
ported measurements of the yield of individual grains of 
silica, olivine, and graphite with 0.09-5 yum radii for 8- 
10 eV photons. Their measurements indeed show larger 
yields than those of the bulk materials. However, the mea- 
surements do not agree with the yield enhancement fac- 
tors adopted in WD01 and W+06 models accounting for 
the small-size particle effect qualitatively as well as quan- 
titatively. Clearly, we need more experiments and theo- 
retical investigations of the photoelectric yield in future. 

2.1.3 Equilibrium charge 

We need to specify the radiation filed incident on grains 
in the IGM: the cosmic background radiation. We assume 
a simple prescription of the radiation. The intensity of 
the radiation at the Lyman limit is estimated from obser- 
vations of the proximity effect and the Lyman a forest 
opacity (e.g., Scott et al. 2000). A typical value of the 
intensity at the Lyman limit is Jl = 1 x 10~ 21 erg s _1 
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cm -2 Hz -1 sr -1 (e.g., Scott et al. 2000). We simply as- 
sume a power-law as the spectral shape: J v oc v~ v . A 
typical value of p is unity (e.g., Haardt & Madau 1996). 
With such a radiation filed, the grains in the IGM are pos- 
itively charged. 

A typical charging time-scale is very short. For ex- 
ample, the collisional charging rate of electron, R c ~ 
5.6 x 10~ 6 s- 1 for n c = 10~ 5 cm- 3 , T = 10 4 K, a = 0.1 
pm, and Z d = 1700, which is the equilibrium charge 
of graphite or silicate grains for these parameters and 
J v = IQ- 21 ^/^)- 1 erg s" 1 cm~ 2 Hz T 1 sr" 1 . Thus, 
the typical charging time-scale ist~l/i? c ^6x 10~ 3 
yr. Therefore, the grain charge can be in equilibrium. We 
set dZ d /dt = in equation (1) and obtain the equilibrium 
charge of the IGM grains. 

2.2 Heating rates 

2.2.1 Heating rate per a grain 

The net heating rate per a grain with the radius a is ex- 
pressed as (e.g., Weingartner & Draine 2001) 

7(a) = i?pe£pe(a)-|i? e |£eCr), (6) 

where £ pc (a) is the mean kinetic energy of photoelectrons 
from a grain with the radius a and £ e (T) is that of elec- 
trons colliding with the grain. The second term accounts 
for the cooling by the electron capture. If we assume 
the Maxwellian velocity distribution for the electrons, 
£ e {T) = k B T(2 + <j>)/{\ + 4>), where = Z d e 2 /ak B T 
(for Z d > 0; Draine 1978). We note that £ c is - 1% of 
£ po in the current setting. 

The mean energy of the photoelectrons is given by 

7ra 2 /' 1/max 47r 7 

£M = -5~ Q„(a)YE„(a,Z d )—^dv, (7) 

H pe J IIV 

and 

YE v {a, Z d ) = ]T Y u k (a, Z d ){E c ) k v {a, Z d ) , (8) 
k 

where Y k is the photoelectric yield of fc-th emission 
process, e.g., primary electrons from the band structure, 
Auger electron, etc., and {E e ) k is the mean energy of 
electrons emitted by fc-th process with the absorbed pho- 
ton energy hv. The estimation of {E e ) k is based on the 



assumed energy distribution of the electrons. Following 
Weingartner et al. (2006), we adopt a parabolic function 
for the primary and the auger electrons and a function in- 
troduced by Draine & Salpeter (1979) for the secondary 
electrons, which were derived to fit some experimental re- 
sults. 

In the IGM, the grains are positively charged. Then, 
the proton collisional charging rate is negligible. Thus, 
the photoelectric charging rate balances with the electron 
collisional charging rate: R pc + R c = 0. In this case, 
equation (6) is reduced to 

7(a) = |i? e |(£pe - £e) 

2 (8k B T\ 1/2 / eV d \ 

* ^H^rJ ter-' (9) 

where we have used equations (2-A) for R e and Vd = 
Z d e/a is the grain electric potential {eV d /k B T 3> 1 and 
£p C » £ c for the IGM). As found later in Fig. 2, the 
electric potential depends on the grain size weakly in the 
W+06 yield case. We confirmed that the mean energy of 
photoelectrons also depends on the grain size weakly. As 
a result, the heating rate per a grain is roughly propor- 
tional to the square of the size, which is shown later in 
Fig. 3. 

2.2.2 Total photoelectric heating rate 

To estimate the total photoelectric heating rate per unit 
volume, we need to specify the amount and the size dis- 
tribution of dust grains. A power law type distribution 
for grain size is familiar in the interstellar medium of the 
Milky Way since classical work by Mathis, Rumpl, & 
Nordsieck (1977; hereafter MRN). Power law is expected 
to be achieved by resulting from coagulation, shattering, 
and sputtering processes (e.g., Jones, Tielens, & Hollen- 
bach 1996). Here we express the power law distribution 
as n(a) = Aa~ q , where n(a)da is the number density 
of grains with the radius between a and a + da. For the 
MRN distribution, q = 3.5 (see Table 1). The normal- 
ization A is determined from the total dust mass density 
p d = J" max m{a)n{a)da, where m{a) = (47r/3)g»a 3 is 
the mass of grains with the radius a, g (~ 3 gem -3 ) is 
the grain material density, a m j n and a max are the mini- 
mum and maximum radius, respectively. The dust mass 
density pd is given by p d = m v nB\D, where m p is the 
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proton mass, «h is the hydrogen number density, and V is 
the dust-to-gas mass ratio. We assume V = 10~ 4 , which 
is about two orders of magnitude smaller than that in the 
Milky Way's ISM. Then, the total photoelectric heating 
rate is 



pe 



j(a)n(a)da . 



(10) 



Let us consider a typical size for the total heating rate. 
Using the grain number density «d = J a Qmax n(a)da, 
we can define a mean heating rate per a grain as (7) = 
la l( a ) n ( a )da/rid and a mean mass per a grain as 
(ma) = Pd/n-d = Xf max m(a)n{a)da/nd- Then, 
the total heating rate is reduced to r po = (7) rid = 
(j) Pd / (nid) ■ The heating rate per a grain can be approx- 
imated to 7(a) « 70a 2 as seen in §2.2.1 (see also Fig. 3 
and §3.1.2), where 70 is a normalization. The grain mass 
is m(a) = (47r/3)pa 3 . Then, we obtain 



r _ 3p d 7o 
where a typical size (a) is given by 



f amax a 3 n(a)da 
( a \ = Ja ^in l_L 



(11) 



(12) 



Note that a larger typical size results in a smaller total 
heating rate because of a smaller number density of grains 
for a fixed dust mass. 

2.2.3 Photoionization heating rates 

For comparison with the photoelectric heating rate by 
grains, we estimate the photoionization heating rates of 
hydrogen and helium. The net HI photoionization heating 
rate is 



r pi — n m R pi £ pi - n m iR rc £, 



where R™ — f°£ T o-™A-Kj v /hvdv is the HI photoioniza- 
tion rate, R™ = n c a™(T) is the HI recombination rate, 
£ P T = (V^pi 1 ) Lh, o?4xJ v /hv(hv - hv™)dv is the 
mean kinetic energy of the HI photoionized electrons, <r™ 
is the HI photoionization cross section, i>™ is the HI Ly- 
man limit frequency, uhi, tihii, and n c are the neutral hy- 
drogen, ionized hydrogen, and electron number densities, 



respectively, a™ (T) is the Case A HI recombination co- 
efficient for the gas temperature T (Osterbrock & Ferland 
2006), and £ gas is the mean kinetic energy lost from the 
gas per one recombination. If we assume that J v oc v~ v 
and cr™ oc we obtain £™ = hu^ l /(p + 2). If we 
take into account the gas cooling by free-free emission, 
£ g as ~ k B T for the Case A and T = 10 4 K (Osterbrock 
& Ferland 2006). If we assume the ionization equilib- 
rium, riHi-R™ = riHii-R™, we obtain 



^HI 



(14) 



where we have assumed nnn = n c = hh with rin being 
the hydrogen number density, that is, the neutral fraction 
is assumed to be very small. The net Hell photoionization 
heating rate is likewise 

- nHcnH^ ClI (r)(fp1 Cl1 - £gas), (15) 

n Ho is the helium number density, a^ eI1 (T) is the Hell re- 
combination rate, and ^ cI1 is the mean kinetic energy of 
the Hell photoionized electrons. We assume nne/^H = 
0.1. 



3 Results 

3.1 Comparison between the two yield 
models 

We compare the grain charge and heating rates with the 
WD01 and W+06 models quantitatively in the IGM en- 
vironment. Weingartner et al. (2006) showed the grain 
charges in QSO environments which are similar situation 
with a similar radiation field in this paper. However, they 
did not show the heating rates in the environment. 



(13) 3.1.1 Electric potential 



In Fig. 2, we compare the electric potentials of the W+06 
model (solid lines) with those of the WD01 model (dashed 
lines). We showed two cases of the spectrum of the radia- 
tion field; one has a hard spectrum as a background radia- 
tion dominated by QSOs which is the case with the spec- 
tral index p = 1, and the other has a soft spectrum with 
p = 5 for a comparison. Other assumed quantities are ap- 
propriate for the IGM at the redshift z ~ 3 and shown in 
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the panels. The radiation fields assumed here correspond 
to the ionization parameter U = n,i on /nu, which is the 
number density ratio of ionizing photons and hydrogen 
nucleus, of 6.3 for p = 1 and of 1.3 for p = 5. Weingart- 
ner et al. (2006) showed electric potentials in their Figs. 6 
and 7 with U = 0.1-100. We find that our calculations 
are quantitatively well matched with theirs. 

We find in Fig. 2 that for the hard spectrum case, the 
grain electric potentials with the W+06 yield model are 
much smaller than those with the WD01 model, especially 
for larger grain sizes. On the other hand, for the soft spec- 
trum case, the difference is very small, less than 4%. This 
is because the main difference between the W+06 yield 
and the WD01 yield is found in the primary photoelec- 
tron yield at <~ 100 eV due to the photon/electron transfer 
in a grain as shown in Fig. 1 . In the soft spectrum case, 
since there are not many photons around the energy, we 
do not find a significant difference between the two yield 
models. For smaller grain sizes, the yield reduction by 
the photon/electron transfer is small as found in Figs. 4 
and 5 of Weingartner et al. (2006). Thus, we do not find a 
significant difference in the electric potentials for smaller 
grain sizes in Fig. 2, either. 

The electrostatic stress on a grain may cause the grain 
destruction by the Coulomb explosion (e.g., Draine & 
Salpeter 1979). The critical electric potential is V ma x = 
1063 V (Sd/10 10 dyn cm- 2 ) 1 / 2 ( a /o.l M m), where S d 
is the tensile strength of grains, which is very uncertain. 
Perfect crystal structure may have Sd ~ 10 11 dyn cm -2 
(Draine & Salpeter 1979), but imperfections would re- 
duce the strength as Sd ~ 10 10 dyn cm -2 (Fruchter et 
al. 2001). Following Weingartner et al. (2006), we show 
two cases of the critical potential with Sd ~ 10 10 and 
10 11 dyn cm -2 in Fig. 2 as the dotted lines. The critical 
potential by the ion field emission is similar to the case 
with S d ~ 10 11 dyn cm~ 2 (Draine & Salpeter 1979). We 
find that grains smaller than 20-30 A in the hard radiation 
field may be destroyed by the Coulomb explosion. Then, 
there may be no very small grains in the IGM. 

3.1.2 Photoelectric heating rate 

Fig. 3 shows the photoelectric heating rate per a grain in 
a typical z ~ 3 IGM environment with a hard radiation; 
graphite grains are shown by solid lines and silicate grains 
are shown by the dashed lines. In the panel (a), we show 



the absolute value of the heating rate for the W+06 yield 
model. As expected in equation (9), the heating rate is 
nicely proportional to a 2 , square of the size. However, 
the slope becomes gradually steep for a small (< 100 A) 
grain size. 

In the panel (b), we compare the two heating rates with 
the W+06 model and the WD01 model. We find that the 
heating rate with the W+06 model is much smaller than 
that with the WD01 model: a factor of 10 smaller for the 
largest grain size. This is because (1) reduction of the 
grain electric potential and (2) reduction of the mean en- 
ergy of the photoelectron in the W+06 model. As found 
in equation (9), the heating rate per a grain is proportional 
to the product of the potential and the mean photoelec- 
tron energy. As seen in Fig. 2, the W+06 model expects 
smaller potential because of the reduction of the yield at 
<~ 100 eV. The yield reduction also causes the reduction 
of the mean energy of the photoelectron as expected in 
equation (8). Therefore, we have up to about a factor of 
10 reduction of the heating rate with the W+06 model. 

Fig. 4 shows a comparison of the total heating rates by 
the W+06 model (solid lines) and by the WD01 model 
(dotted lines). The horizontal axis is the assumed hydro- 
gen number density. We also show the redshift at which 
the number density on the horizontal axis corresponds to 
the mean density of the Universe. We have assumed the 
MRN grain size distribution (see Table 1). We find that 
the total heating rate with the W+06 yield is a factor of 
2-4 smaller than that with the WD01 yield. 

For a comparison, we also show the HI and Hell pho- 
toionization heating rates in Fig. 4. We have assumed the 
ionization equilibrium for them. When we assume the 
dust-to-gas ratio in the IGM is 1 % of that in the Milky 
Way (i.e. V = 10~ 4 ), the dust photoelectric heating 
dominates the HI and Hell photoionization heatings if the 
hydrogen number density is less than 10~ 6 -10~ 5 cm~ 3 
which correspond to the redshift z ~ 1-2. We note that 
the dust heating is the most important mechanism in the 
IGM at z = even with the W+06 yield model if the IGM 
has dust with 1 % dust-to-gas ratio of the Milky Way and 
with the MRN size distribution. 

3.2 Effect of the grain size distribution 

The size distribution of the intergalactic dust grains 
should be important for the photoelectric heating rate via 
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the typical size defined by equation (12). However, it is 
quite uncertain. Thus, we examine several possibilities of 
the size distribution in this section. Table 1 is a summary 
of the size distribution considered here. 

The grain size distribution in the Milky Way has been 
approximated to be a power-law since Mathis et al. (1977) 
suggested as n(a) oc a~ q with q = 3.5. This MRN distri- 
bution is a reference case and is already adopted in Fig. 4. 
During the grain transport from galaxies to the IGM, there 
may be size filtering mechanisms. For example, Ferrara 
et al. (1991) showed that sputtering in the hot gas filling 
the galactic halo efficiently destroys grains smaller than 
^0.1 /im. Bianchi & Ferrara (2005) also showed that 
only grains larger than ^0.1 ^m reach a significant dis- 
tance (a few x 100 kpc) from the parent galaxies by calcu- 
lating the grain ejection by the radiation pressure and the 
grain destruction by the sputtering simultaneously. Here, 
we consider a simple size distribution of the MRN with 
> 0.1 /im grains as the BF05 model. 

In the early Universe, the dominant source of dust 
grains is different from that in the current Milky Way. Al- 
though asymptotic giant branch stars are considered to be 
the main dust source in the Milky Way (e.g., Dwek 1998), 
there is not enough time for stars to evolve to the phase 
in the early Universe at the redshift z > 6. However, a 
plenty of dust is found in QSOs at z > 6 (Bertoldi et 
al. 2003). SNe are the candidate of the dust source in the 
early Universe (e.g., Nozawa et al. 2003) and the observed 
extinction curve of dust associated with the QSO is com- 
patible with those expected from the grains produced by 
SNe (Maiolino et al. 2004, Hirashita et al. 2005,2008). 
Thus, we consider the size distribution expected from the 
SNe dust production model by Nozawa et al. (2003) as 
the N03 model. In addition, we adopt the size distribution 
expected by Nozawa et al. (2007), who explored the ef- 
fect of the dust destruction by the reverse shock in the SN 
remnant, as the N07 model. 

Finally, we adopt a hypothetical size distribution con- 
sisting of only small grains as a comparison case; the 
MRN distribution with the maximum size of 250 A as the 
SG (small grain) model. 

Fig. 5 shows a comparison of total heating rates with 
the five size distributions considered here. All the cases 
are assumed the W+06 yield model and physical con- 
ditions appropriate for the IGM. The case of the BF05 
model (triple-dot-dashed line) is a factor of ~ 5 smaller 



than that of the MRN model (thick solid line). This re- 
duction factor is simply accounted for by the ratio of the 
typical sizes of the two models: 0.16 [im for the BF05 
model and 350 A for the MRN model (see Table 1). The 
same thing is true for the N07 model (dotted line) and the 
SG model (thin solid line). The result of the N03 model 
(dashed line) coincides with that of the MRN model be- 
cause their typical sizes are similar. In any case, we have 
a smaller number of grains for a larger typical size if the 
total dust mass is fixed. Then, the heating rate is reduced. 
We note that the dust photoelectric heating is still domi- 
nant or important mechanism relative to the HI and Hell 
photoionization heatings in the z = IGM even with the 
BF05 model if the dust-to-gas ratio in the IGM is 1 % of 
that in the Milky Way. 

3.3 A simple formula of the dust photoelec- 
tric heating rate 

Fig. 6 shows the effect of different settings of the cal- 
culation on the dust photoelectric heating rate: (a) vari- 
ous intensities of the background radiation and (b) vari- 
ous temperatures of the gas. The W+06 yield model and 
the MRN size distribution are assumed. We also assume 
that the dust consists of a mixture of graphite and sili- 
cate with the mass ratio of 1 : 1 . The spectral index of the 
background radiation is always set unity. In the weak- 
est intensity case (squares in the panel [a]), the equilib- 
rium charges for smallest grains (< 0.01 /im) are less 
than 3 in the electric charge unit for the hydrogen den- 
sity jih > 2 x 10 -5 cm~ 3 . In these cases, the effect of 
the "image potential" (Draine & Sutin 1987) is not neg- 
ligible, then, the current calculations are no longer valid. 
We note that all the cases shown in Fig. 6 have an equilib- 
rium charge much larger than 3 for all grains in the size 
distribution. 

The resultant heating rates are well expressed as 



pe 



1.2 x 10~ 34 ergs" 1 cm" 3 



V 



( nH V 
V10~ 5 cm- 3 / 



4/3 / t 



10 4 K 



10- V \10- 5 cm- 3 

Jl 

1Q- 21 erg s _1 cm- 2 Hz- 1 sr- 1 



-1/6 



2/3 



(16) 



which is shown in Fig. 6 as solid lines. The indices 
in this formula can be derived analytically follow- 
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ing Inoue & Kamaya (2004). From eqs. (A4) and 
(A7) in Inoue & Kamaya (2004), we find r pc oc 

j2/(p+/3+l) w 2-2/(p+/3+l) T 3/2_(2p+2/3+l)/(p+/3+l) 
L H ' 

where (3 is the emissivity (or absorption) index of the 
dust: Q v oc v~P . Here, we have p = 1 and (3 « 1, then, 
we obtain the indices in eq. (16). 

The deviation of the heating rates from the formula for 
T = 10 5 K and n H > 2 x 10~ 5 cm~ 3 is due to the 
relative significance of the cooling by the electron capture 
(see eq. [6]). Indeed, we find that the mean energy of 
photoelectrons from small (< 0.01 fim) graphite grains 
is smaller than the mean kinetic energy of the 10 5 K gas 
in the case of iih = 10~ 4 cm~ 3 , then, the net heating 
rate per such a graphite grain is negative. Thus, we have 
a reduction of the total heating rate for T = 10 5 K and 
tin > 2 x 10~ 5 cm~ 3 as found in Fig. 6 (b) although the 
heating rate by silicate grains is still positive. 

The validity of the formula presented in equation (16) 
is ensured for n H = 10~ 7 -10~ 4 cm~ 3 , J L = 10~ 23 - 
10~ 20 erg s- 1 cm" 2 Hz 1 sr"\ and T = 10 3 -10 5 K 
within a uncertainty of 30%, except for n H > 2 x 10~ 5 
cm -3 with Jl = 10~ 23 erg s _1 cm~ 2 Hz -1 sr _1 or T = 
10 5 K. Note that there may be much larger uncertainty in 
the photoelectric yield model. If one likes another size 
distribution rather than the standard MRN, for example 
the BF05 model discussed in §3.3, the heating rate might 
be scaled by a factor found in Fig. 5 or the ratio of the 
typical sizes in Table 1 . 

4 Discussions 

4.1 Amount of the intergalactic dust 

Inoue & Kamaya (2003,2004) discussed the effect of the 
photoelectric heating by the intergalactic dust on the ther- 
mal history of the IGM, and then, obtained an upper limit 
of the intergalactic dust amount. However, we have al- 
ready seen that the W+06 yield model results in a factor 
of 2-4 reduction of the photoelectric heating rate relative 
to the WD01 model which was adopted in Inoue & Ka- 
maya (2003,2004). We can conclude that the upper limits 
obtained from the IGM thermal history are raised by a 
few factor. Even in this case, the final limit obtained by 
Inoue & Kamaya (2004), which is that the intergalactic 
dust mass should be less than 10% of the metal mass pro- 



duced in galaxies, is not affected because the limit was ob- 
tained mainly from the reddening measurements of SNe 
la at z — 0.5, especially for ~ 0.1 ^m size grains. 

4.2 Can grains cause an inverted 
temperature-density relation in the 
IGM? 

Bolton et al. (2008) recently suggest an inverted 
temperature-density relation in the low density IGM at 
z = 2-3. Previously the temperature in the low density 
IGM is thought to be proportional to the density positively 
(e.g., Hui & Gnedin 1997). However, Bolton et al. (2008) 
examined carefully the probability distribution function 
(PDF) of the flux in QSOs' spectra through the Lyman a 
forest in the IGM, and found that the observed PDF is ex- 
plained better by the negatively proportional temperature- 
density relation; lower density IGM is hotter. This needs 
a more efficient heating source for lower density IGM. 
Bolton et al. (2008) suggested a radiation transfer effect 
(e.g., Abel & Haehnelt 1999) for the mechanism. 

The intergalactic dust may contribute to the heating in 
the low density IGM. As shown in Figs. 4 and 5, the 
importance of the dust photoelectric heating increases in 
lower density medium which is plausible for the inverted 
temperature-density relation. For example, we expect a 
factor of ~ 2 larger heating rate by dust than Hell pho- 
toionization heating in a medium with 1/10 of the mean 
density at z — 2 for the MRN size distribution and 1% 
dust-to-gas ratio of the Milky Way. Thus, the dust pho- 
toelectric heating may cause the inverted temperature- 
density relation observed in the low density IGM at z = 
2-3. This point should be examined further by imple- 
menting the dust heating in a cosmological hydrodynam- 
ics simulation. For this, the formula presented in equa- 
tion (16) will be useful. 

4.3 Photoelectric effect before the cosmic 
reionization 

Finally, we examine if the dust photoelectric heating is 
efficient in the IGM before the cosmic reionization. Be- 
cause of the prominent Gunn-Peterson trough in QSOs' 
spectra (e.g., Fan et al. 2006), the cosmic reionization 
epoch should be at z > 6. Here we consider the IGM 
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at z ~ 10. 

Before the reionization, the ionizing background ra- 
diation does not exist although a nonionizing UV back- 
ground can be established by primordial galaxies or active 
blackhole-accretion disk systems. An X-ray background 
radiation may also exist (e.g., Venkatesan et al. 2001). We 
consider two cases; one is the case with only nonionizing 
UV background radiation and the other is the case with 
additional X-ray background radiation. For simplicity, 
we assume the background radiation to be a power-law 
with the spectral index p = 1 and the intensity at the Ly- 
man limit J L = 1 x 10~ 21 erg s _1 cm~ 2 Hz 1 sr _1 . 
However, we assume no intensity between -EJjy* ~ 13.6 
eV and Ef in = 300 eV. Thus, in the nonionizing UV 
only case, we have the background radiation only below 
£"uy x = 13.6 eV. In the case with the X-ray background, 
we have radiation below E™y x = 13.6 eV and above 
E™ in = 300 eV. The dust-to-gas ratio in the IGM at 
z ~ 10 is of course unknown, but we assume 1% dust-to- 
gas ratio of the Milky Way as an example, i.e. V = 10~ 4 . 
Note that the results obtained in the following discussions 
are linearly scaled by the value of V. The mean hydrogen 
density in the Universe at z <~ 10 is 3 x 10~ 4 cm~ 3 . Ta- 
ble 2 is a summary of the assumed quantities and results 
obtained below. 

In the nonionizing UV radiation only case, there is no 
efficient heating mechanism for the whole of the Universe 
although primordial objects can heat up their surrounding 
gas locally. Thus, the temperature of the gas far away 
the sources is kept to be that of the cosmic background 
radiation at the epoch: ~ 30 K. The electron fraction 
x c , i.e. the number density of electron relative to that 
of hydrogen nucleus, is ~ 10~ 4 in this low temperature 
IGM (Galli & Palla 1998). The nonionizing UV photons 
still cause the photoelectric effect of grains. In the as- 
sumed setting, we have found that grains are positively 
charged and the dust photoelectric heating rate becomes 
r po ~ 7 x 10~ 36 erg s _1 cm -3 fortheMRN size distri- 
bution with a graphite and silicate mixture (50% each in 
mass). We compare this heating rate with the gas thermal 
energy density: U gas — (3/2)nnkBT. The time-scale 
doubling the gas temperature with the photoelectric heat- 
ing is given by t pc = U gas /T pc ~ 9 x 10 9 yr Since the age 
of the Universe at z = 10 is about 5 x 10 8 yr, we conclude 
that the dust photoelectric heating is not very efficient in 
this case although it may be the strongest heating mecha- 



nism for the IGM. 

In the case with the additional X-ray background ra- 
diation, the IGM is partially ionized by the X-ray and 
the temperature becomes ~ 10 4 K (e.g., Venkatesan et 
al. 2001). If we assume the ionization equilibrium and 
optically thin for the X-ray, the electron fraction becomes 
x c ~ 0.3 for the current setting of the X-ray background. 
In this medium, the grains are positively charged and the 
dust photoelectric heating rate becomes r po ~ 2 x 10~ 33 
erg s~ 4 cm -3 . We have assumed the MRN size distri- 
bution with a graphite and silicate mixture (50% each in 
mass), again. However, the HI photoionization heating is 
much more efficient as r™ x ~ 2 x 10~ 30 erg s _1 cm~ 3 . 
Therefore, we again conclude that the dust photoelectric 
heating is negligible in the early Universe filled with an 
X-ray background radiation. 

5 Conclusion 

We have updated our calculations made in Inoue & Ka- 
maya (2003,2004) of the dust photoelectric heating in the 
IGM with the new model of the dust photoelectric yield 
by Weingartner et al. (2006). This new yield model takes 
into account the effect of the photon and electron trans- 
fer in a grain, the photoelectric emission from inner shells 
of grain constituent atoms, the Auger electron emission, 
and the secondary electron emission. A comparison with 
the previous yield model by Weingartner & Draine (2001) 
show that the new yield is smaller than the old one for 
~ 100 eV photons. This reduction of the yield is due to 
the photon/electron transfer effect, and reduces the elec- 
tric potential on grains and the heating rate significantly. 
For example, if we integrate over the grain size with the 
standard MRN distribution, the dust photoelectric heating 
rate with the new yield model is a factor of 2-4 smaller 
than that with the old yield model. The photoelectric heat- 
ing rate is more important in lower density medium. If 
the dust-to-gas ratio in the IGM is 1 % of that in the Milky 
Way and the size distribution is the standard MRN model, 
the dust heating rate dominates the HI and Hell photoion- 
ization heating rates when the gas number density is less 
than ~ 10~ 6 cm~ 3 even with the new yield model. 

We have examined the effect of the size distribution 
function on the heating rate because the heating rate is 
inversely proportional to the typical grain size as found in 
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equation (11). Bianchi & Ferrara (2005) suggested that 
the size of the intergalactic dust is larger than ~ 0.1 
because smaller grains are destroyed by sputtering in the 
hot gas halo during the transport of grains from the parent 
galaxy to the IGM. In this case, the heating rate is re- 
duced by a factor of ~ 5 relative to that with the standard 
MRN size distribution. The size distributions expected 
by the dust formation model in supernova ejecta are also 
examined. The heating rate with the size distribution of 
the grains just produced in the ejecta is very similar to 
that with the MRN distribution. In contrast, the heating 
rate with the size distribution of the grains processed by 
the reverse shock in the supernova remnant is a factor of 
~ 3 smaller than that with the MRN model. The shock 
processed grains have larger size than the pristine ones 
because smaller grains are destroyed. On the other hand, 
if we put only small grains in the IGM, the heating rate 
increases significantly. Therefore, we conclude that the 
size distribution of grains in the IGM is an essential pa- 
rameter for determining the dust heating efficiency. Even 
in the worst case considered here, the dust heating is ex- 
pected to be the dominant heating mechanism in the IGM 
at z = if the dust-to-gas ratio in the IGM is 1 % of that 
in the Milky Way. 

Since the dust photoelectric heating rate with the new 
yield model is reduced by a factor of 2-4 relative to that 
with the old yield model, the upper limit on the amount 
of the intergalactic dust obtained by Inoue & Kamaya 
(2003,2004) may be affected. Indeed, the limit based on 
the thermal history of the IGM should be raised by a few 
factor. However, their final upper limit is mainly obtained 
from the reddening measurements of z = 0.5 supernovae 
la. Therefore, their conclusion would not be affected very 
much. 

Bolton et al. (2008) suggested an inverted temperature- 
density relation in the lower density IGM at z = 2-3 
based on recent observations of the Lyman a forest in 
QSOs' spectra. To explain this interesting phenomenon, 
we need a heating mechanism more efficient in a lower 
density medium. The dust photoelectric heating has such 
a property. Indeed, the dust heating rate even with the new 
yield model is a factor of 2 larger than the Hell photoion- 
ization heating rate in a medium with a density of 1/10 of 
the mean in the Universe at z = 2 if the dust-to-gas ratio 
is 1% of that in the Milky Way. Thus, the possibility of the 
dust heating is worth examining more in detail. For this 



aim, the simple formula describing the dust photoelectric 
heating in the IGM presented in equation (16) will be very 
useful. 

Finally, we have discussed the effect of the dust photo- 
electric heating in the early Universe. Before the cosmic 
reionization, the ionizing background radiation is not es- 
tablished, but there may be nonionizing UV background 
and X-ray background radiations. In the low tempera- 
ture IGM only with a nonionizing UV background radi- 
ation, the dust photoelectric heating is not very efficient 
although it may be the strongest heating mechanism in 
the medium. In the partially ionized IGM with an X-ray 
background radiation, the HI photoionization heating rate 
is three orders of magnitude larger than the dust heating 
rate if the dust-to-gas ratio is 1% of that in the Milky Way. 
Therefore, we conclude that the dust photoelectric heating 
in the early Universe is not very important at least in the 
mean density environment. 
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Figure 1: Photoelectric yield models of 0.1 fim neutral 
(a) graphite and (b) silicate grains. The dotted lines are 
the WD01 model (Weingartner & Draine 2001) and the 
solid lines are the W+06 model (Weingartner et al. 2006). 
The W+06 model consists of three processes: the primary 
photoelectron emission (short dashed line), the secondary 
electron emission (long dashed line), and the Auger elec- 
tron emission (dot-dashed line). 
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Figure 2: Equilibrium electric potential as a function of 
grain size: (a) graphite and (b) silicate. The solid lines 
are the W+06 model and the dashed lines are the WD01 
model. The thick lines are the case with the spectral index 
of the radiation field p = 1 and the thin lines are the case 
with p = 5. Other assumed quantities are noted in the 
panels as the hydrogen density tih,-5 = "-h/10~ 5 cm -3 , 
the gas temperature T4 = T/10 4 K, and the radiation 
intensity Jl.-21 = Jl/10~ 21 ergs -1 cm -2 Hz -1 sr _1 . 
The dotted lines show the critical electric potential where 
the grain destruction occurs by the Coulomb explosion; 
the upper lines are the case with the tensile strength of 
10 11 dyn cm -2 and the lower lines are the case with 10 10 
dyn cm -2 . 
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Figure 3: Photoelectric heating rate per a grain as a func- 
tion of grain size: (a) heating rate with the W+06 yield 
model and (b) ratio of the heating rate with the W+06 
yield model to that with the WD01 yield model. The solid 
lines are the graphite case and the dashed lines are the sil- 
icate case. The assumed quantities are noted in the panels 
as the hydrogen density njj.-s = rijj/10 -5 cm -3 , the 
gas temperature T4 = T/10 K, the radiation intensity 
Jh-21 = Jl/10 -21 ergs -1 cm -2 Hz -1 sr _1 , and the 
spectral index of the radiation field p = 1, The thin solid 
line in the panel (a) shows the slope proportional to the 
square of size. 
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Figure 4: Photoelectric heating rates as a function of hy- 
drogen number density: (a) graphite and (b) silicate. The 
solid lines are the W+06 model and the dotted lines are the 
WD01 model. The assumed grain size distribution is the 
so-called MRN distribution (Mathis, Rumpl, & Nordsieck 
1977). Other assumed quantities are noted in the panels as 
the gas temperature T4 = T/10 4 K, the radiation inten- 
sity Jl.-21 = Jl/10 -21 ergs -1 cm -2 Hz -1 sr -1 , the 
spectral index of the radiation field p, and the dust-to-gas 
mass ratio 25,-4 = 2?/10 -4 . The dashed lines are the 
HI photoionization heating rate and the dot-dashed lines 
are the Hell photoionization heating rate with assuming 
the ionization equilibrium. We also show the redshift at 
which the number density on the horizontal axis corre- 
sponds to the mean density of the Universe. 
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Table 1: Possible size distributions of the intergalactic 
dust. 



MRN Mathis, Rumpl, & Nordsieck (1977) 

single power law a 

q 3.5 

a m j n 50 A 

a max 0.25 /jm 

(a) 350 A 

BF05 Bianchi & Ferrara (2005) 

single power law" 

q 3.5 

dmin 0.1 /Lim 

a ma x 0.25 /im 

(a) 0.16 /im 

N03 Nozawa et al. (2003) 

double power law fc 

gi (a < a c ) 2.5 

52 (a > a c ) 3.5 

Qniin 2 A 

Omax 0.3 fim 

a c 0.01 /im 

(a) 290 A 

N07 Nozawa et al. (2007) 

double power law fc 

qi (a < a c ) 1.0 

g 2 (a > a c ) 2.5 

Clmin 10 A 

Omax 0.3 

a c 0.01 /im 

(a) 0.12 

SG — 
single power law a 

q 3.5 

50 A 

a max 0.025 /xm 

(a) 110 A 




10" 7 10" 6 10" 5 10" 4 

Hydrogen number density [cm -3 ] 

Figure 5: Same as Fig. 4 but for various size distribution 
functions with the W+06 yield model: (a) graphite and 
(b) silicate. The thick solid lines are the MRN case. The 
short-dashed lines are the size distribution expected from 
the grain formation model in supernova ejecta by Nozawa 
et al. (2003). The dotted lines are the size distribution 
expected after the grain destruction by the reverse shock 
in the supernova remnant by Nozawa et al. (2007). The 
triple-dot-dashed lines are the MRN but only size larger 
than 0. 1 /im because of a filtering effect in the transfer of 
grains from galaxies to the IGM suggested by Bianchi & 
Ferrara (2005). The thin solid lines are the MRN but only 
size smaller than 250 A as a comparison. The dashed and 
dot-dashed lines are HI and Hell heating rates. 



a The grain size distribution n(a) oc a~ q . 

b The grain size distribution n(a) oc aT qi for a < a c and 

oc a~ 92 for a > a c . 
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Figure 6: Same as Fig. 4 but for various settings. The pho- 
toelectric yield model is the W+06 model. We assume that 
the dust consists of a mixture of graphite and silicate (50% 
each in mass) with the MRN size distribution, (a) Differ- 
ent intensities at the Lyman limit of the radiation field: 
Jl/10 -21 erg s _1 cm" 2 Hz 1 sr _1 = 10 (diamonds), 1 
(circles), 0.1 (triangles), and 0.01 (squares), (b) Differ- 
ent temperatures of the gas: T/10 4 K = 10 (diamonds), 
1 (circles), and 0.1 (triangles). Other assumed quantities 
are noted in the panels. See the caption of Fig. 4 for the 
notation. The solid lines are the simple formula shown in 
eq. (16). 



Table 2: Photoelectric heating in the early Universe. 



Common setting 


z 


10 






«H 


3 x 10" 4 crrr 


-3 




V 


io- 4 






size distribution 


MRN 






Jl 1 x 


10 21 erg s 1 cm 2 Hz 


- 1 sr- 1 


V 


1 






rp max 


13.6 eV 






Nonionizing UV only 


T 


30 K 






Xq 


10~ 4 








7 x 10~ 36 erg s" 1 


cm 


-3 




9 x 10 9 yr 






With X-ray background 




300 eV 






T 


10 4 K 






X e 


0.3 






r 

x pe 


2 x 10" 33 erg s" 1 


cm 


-3 


pHI 

pi,X 


2 x 10~ 30 erg s" 1 


cm" 


-3 



16 



